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We consider a modified "Swiss cheese" model in Brans-Dicke theory, and use it to discuss the 
evolution of black holes in an expanding universe. We define the black hole radius by the Misner- 
Sharp mass and find their exact time evolutions for dust and vacuum universes of all curvatures. 
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CN ■ I. INTRODUCTION 

1^ ' The evolution of primordial black holes in scalar-tensor theories has been studied in the literature by several authors 
I0~§| ■ Some questions about what happens to a black hole in an expanding universe when G varies in spacetime in these 
00 theories were posed by Barrow [0, who considered two possible scenarios: (a) the effective gravitational "constant", 
G{t), at the black hole horizon changes along with its cosmological evolution so that the size of a black hole is 
approximated by i? = 2G{t)M. (b) G remains constant at the black hole event horizon while it evolves on larger 
' scales; a large inhomogeneity in G is therefore generated. The case (b) was called "gravitational memory" because 
^ , the black hole remembers the value of G at its formation time. In either case the observational constraints on the 
^1 abundance of exploding primordial black holes deduced from the total radiation backgrounds today would be modified 

PsJ ' Scheel, Shapiro and Teukolsky made numerical analyses of dust collapse in Brans-Dicke (BD) theory, showing 
I that the surface area of the event horizon decreases with time, contrary to the case in Einstein theory. Kang H| 
• gave an analytic explanation for the surface-area decrease in BD theory. Later, Jacobson claimed that there is no 
[ "gravitational memory" effect, by analyzing the evolution of a scalar field (j){t, r) in Schwarzschild background He 
O . found a particular solution of the scalar wave equation which matches smoothly between the black hole and a special 
O^' cosmological background and showed that at the event horizon evolves along with its asymptotic value (j){t, oo). 
. However, this solution requires a particular cosmological variation of G{t) to occur in the background and may be 
special. It was also argued that even if the black hole mass in the Einstein frame is constant then its mass in the 
^ ' Jordan frame is time-dependent. 
• '~j In order to investigate the time-dependence of the black-hole mass, Saida and Soda constructed a "cell lattice" 

rN [ universe in BD theory In their model the universe is first tessellated by identical polyhedrons, which are then 
^ ,. replaced by Schwarzschild black holes. It was shown that the black-hole mass has an adiabatic time dependence, 
■ - - ' which is qualitatively different according to the sign of the curvature of the background universe. 

As an extension of Saida and Soda's work, we consider a "Swiss cheese" (or Einstcin-Straus) model in BD theory 
and discuss the evolution of the radius and mass of black holes in an expanding isotropic background universe. The 
usual "Swiss cheese" model refers to a cosmological model in which spherical regions in the Friedmann-Robertson- 
Walker (FRW) universe are replaced by Schwarzschild spacetimes. Here we construct such a model in BD cosmology. 

II. BACKGROUND UNIVERSE IN BRANS-DICKE THEORY 

BD theory is described by the action, 



S = d x\/ —g 



ijj 



(2.1) 
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where </> is BD field, lo is BD parameter, and Cm is the matter Lagrangian. The variations of equation (2.1) with 
respect to g^u and 4> yield the field equations: 



1 7? _ ^'^T ^ 



-TrT. 



2w + 3 

As a background universe, we assume the FRW spacetime: 



ds'^ = -dr + a{ty 



+ r^(dr +sin^ edifi^) } , 



(2.2) 
(2.3) 

(2.4) 



1 — kr^ 

where k = 0,±1 determines the spatial curvature. As an energy-momentum tensor, we introduce a dust fluid: 

Tpi/ = pu^Uu, (2.5) 



where p and are the density and the four- velocity of dust, respectively. The field equations ( p. 2D and (2.3) reduce 
to the following equations for the background universe: 



Sttp 
"30" 



3i/0: 



Snp 



2tj + 3' 

where an overdot denotes d/dt and H = d/a is the Hubble parameter. 

In the flat space case, with fc = 0, we know the general analytic solution for a dust universe ||]: 

a{t) = ao{t - {t ~ i_)^- , <f>{t) ^Mt- t+r+ [t ~ t-r- , 

with constants X± and k± defined by 



A. 



Lo + l±Jl + 2uj/3 l±3Vl + 2w/3 

-, K± = ■ 



3uj + A 



3w-h4 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



and ao, ^io ^^nd t± are arbitrary constants. 

If we take i+ = t-, the general solution ( |2.8| ) reduces to the special power-law solution, (used for example by Saida 
and Soda M): 



a{t) ^ao{t- to) , (j){t) = (jjoit - to) 3-+4 . 



(2.10) 



where to may be set to zero. In the limit of t oo, the general solution (2.8) converges to the special power-law 
solution (2.10). If the present cosmic age t is large enough, observations cannot constrain the relation between t^ and 
t_. Therefore, we keep t+ ~ t^ a free parameter. 

In vacuum case (p — 0) there are analytic solutions for all k. The vacuum solution for fc = is expressed as 



,t— 



with 



1 ± ^1 + 2tj/3 



(2.11) 



(2.12) 



where we have omitted the arbitrary constant to by fixing the origin of the time coordinate t. Introducing the 
conformal time Tj = J dt/a, the vacuum solutions for k = ±1 are expressed as 



k = +1 : a{ri) — (sin?]) 2 (cost]) 2 ^ (^(7^) — (tan?])'*', 
k — —1 : a{r/) = (sinh?7)~2-(cosh77)~5-^ (j)(r]) = (tanh?])'''. 



with 



A = ±- 



'3 + 2uj 



(2.13) 
(2.14) 

(2.15) 
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III. MODIFIED SWISS-CHEESE MODEL 



Now we consider a model for a black hole embedded in the FRW universe. We replace a sphere in the FRW 
spacetime with a vacuum region which contains a black hole. Here "vacuum" means T^^ — 0, and does not imply 
that 7?.^y = due to the existence of BD field. 

Extending Israel's junction conditions for a singular (or regular) hypersurface Sakai and Maeda have studied 
bubble dynamics in the inflationary universe [pl| . It was found that one can solve the equations of motion for the 
boundary without knowing the interior metric if the interior is vacuum, T^^ = 0, or has only vacuum energy (a 
cosmological constant), T^^i, = —pg^^^. Applying this method to the present model, we can determine the mass and 
the radius of a black hole without specifying the interior metric, as we shall show below. 

Let us consider a spherical hypersurface E which divides a spacetime into two regions, (outside) and (inside). 
We define a unit space-like vector, N^j^, which is orthogonal to S and points from V~ to y+. In order to describe 
the behaviour of the boundary, we introduce a Gaussian normal coordinate system, {n,x^) = {n,T,6,ip), where t is 
chosen to be the proper time on the boundary. Hereafter, we denote by the value of any field variable \1/ defined 
on E by taking limits from . 

For the matter field, we consider a dust (or vacuum) medium for V'^ and vacuum for V^: 



0. 



(3.1) 



Although we assume a smooth boundary at which there is no surface density, it is not obvious that this matching is 
possible at all times. Therefore, we introduce a surface energy-momentum tensor on the boundary surface. 



Sij = lim J dn — diag(— cr, tn, w), 



(3.2) 



where a and w denote the surface energy-density and the surface pressure of S, respectively. 

If we introduce the extrinsic curvature tensor of the world hypersurface E, Kij = Ni-j, we can write the junction 
conditions on E as 



47r 



2w + 3 



TrS'7i 



(3.3) 



(3.4) 



2tt ( 



(3.5) 



where we have defined the jump in any quantity ^' by the bracket [4*]^ = 'f^ — ^' and the three-dimensional covariant 
derivative by the vertical bar |. The junction condition for BD field is derived from equation (p. 31) as 



247r 
'3 + 2uj 



TiS, 



(3.6) 



which implies that (f) is continuous at E and inhomogeneous in 

The extrinsic curvature tensor of E in the homogeneous region is given by 11] 



j{l + vHR) _ e 
R 




i?2 i/2 



(3.7) 
(3.8) 



where 



R = a{t)r\^, V = a 



57 



7T=f' ^ ^ sign(i^.^) = sign 



(3.9) 



From equations (3.1), (3.2), (|3.4[), (|3.5D, (3.7)-( |3.9|), we obtain the equations of motion 
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(3.10) 



3 '^^ 



„ , 2w dr 

l-2w]vH - — — 
R dx 



2Txa 



1 + 4w - 



(1 - 2wf \ -i'^v'^p 
{2uj + 3) 



(3.11) 



da 



2(j{l + w) dR 
'dt 



R 



(3.12) 



where w = tu/a. 

Once initial values of R, v, and cr are given, the equations of motion (3. IC)-(O^) determine their evolution. As 
discussed in [0, initial values should satisfy the angular component of (3.3), 



j{l + vHR) 




SiraR f LU + l + w 
2uj + 3 



where we have chosen e = +1. Rms is defined by 

i?Ms = i?"(i-5^"i?-i?:.), 

where R~ is defined as R^ = ^/gee at E on the F^side. Since the Misner-Sharp mass is defined as 

Mms ^ §(i - g'^RT.R:.) - 



(3.13) 



(3.14) 



(3.15) 



we call Rms the "Misner-Sharp radius" . Note that Rms is a purely geometrical quantity and independent of theories 
of gravitation. 

If we considered a spherical bubble in which there is no black hole (or singularit y), w e would have to solve the 
field equations with a regularity condition at the centre and the boundary condition (3.13), as done in However, 
because we are interested in blac k hole solutions, we do not have to take a central regularity condition into account. 
Thus, we can use equation ( ^.13 ) to determine Rms- Of course, our treatment does not guarantee the existence of a 
black hole in the center. We suppose a black hole as the most interesting object, which can be modeled by appropriate 
choices of matter and BD field configuration. 

At the initial time, we suppose v = and a = zu = 0, so Rms is given by 



Rms — R [ H + 



(3.16) 



Let us now discuss whether v and a remain zero during the ensuing evolution. Suppose w — 0, then the only nontrivial 
term in equation (3.11) is ^'^v^pja. If v and a ev olved from zero, then equation (3.13) shows vH ~ o'/'/', so that 
pv^ ja ~ pv/H(j). Therefore, equations ( |3.11 ) and ( 3.12 ) guarantee that, if w = cr = at a certain time, v — a ~ Q 
at all tim e. In terestingly, this result is true only for the dust case, -uj/a = 0; otherwise the term {2w I R){dr / dx) in 
equation ( ^.11 ) would shift v from zero. 

In the case of Schwarzschild spacetime, the Misner-Sharp radius coincides with the event horizon. Although this 
is not necessarily true for general spacetimes, we speculate that the Misner-Sharp radius is a well-defined measure of 
the size of a black hole. In the next section, we calculate the evolution of Rms for black holes in several background 
cosmological models. 



IV. EVOLUTION OF BLACK HOLES 



The evolution of the Misner-Sharp radius for the k — Q dust universe is given by equations (2.8) and (3.16) 



Rms = alrl 



t - U 



t - 1- 



3X- 



(4.1) 



where rg is the comoving radius of the vacuum region. Equation (4.1) shows that the black hole size decreases with 
time. 
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If we define the black hole mass by 



Mms = — ^ — , (4.2) 
then it coincides with the mass defined by Saida and Soda 1^. For the fc = dust universe, we obtain 

mms = ^ + - t+r^+^Ht - t.r-+^-. (4.3) 



It is easy to see that equation (4.3) reduces to Mms = constant, if we choose = which is the same result 
as that found by Saida and Soda ||6[. They also showed Mms increases for fc = +1 and decreases for fc = —1, and 
concluded that the evolution of the mass depends qualitatively on the sign of the curvature of the universe. We should 
note, however, that their conclusion is true only for the special case i+ =t-, or equivalently, only for the asymptotic 
behavior of Mms at i ^ oo. Our results give the general solution for all times. 

Next, let us consider the scalar-field dominated (vacuum) universe. The time-dependent solutions for Rms E^nd 
Mms in Af't, open, and closed universes are given by 

^3 

A: = +1 : Rms = (cos 77) (sin 77)"" (cos 27? - sin 277 - A), (4.6) 

7^ 3-|-A \ I A 

Mms = ^(tan?7)-*'(cos77) ~ {sinvi) 5-(cos 277 - sin 277 - A), (4.7) 

f 3-|-A 1 -|- A 

fc = -1 : Rms = y (coshr;) ^-(sinhr;) 5-(cosh2r7 - sinh27; - A), (4.8) 
Mms = ^(tanh77)^(cosh77) ^-(sinhTy) ^ (cosh 2?; - sinh277 - A). (4.9) 

If we take uj > 500, as applies to the universe today, then both a ^ 0{uj^^) and A ^ 0{uj^^) will be negligible. We 
see that Rms and Mms both decrease with increasing time, except in the contracting phase of the fc = -1-1 universe. 
In a more general scalar-tensor theory, with non-constant w(0) we would expect similar effects to arise and it would 
be possible for large changes in Lu{(f)) to occur in the very early universe despite very slow evolution towards a very 
large value of uj > 1000 today. However, such a calculation should be performed for a black hole in radiation and 
dust-dominated universes. 



V. DISCUSSION 



We have constructed a modified "Swiss cheese" model in Brans-Dicke theory and discussed the evolution of black 
holes for dust and vacuum universes. We defined the size of a black hole Rms by the Misner-Sharp mass, and found 
that this size always decreases as long as the universe is in an expanding phase. Although we have not specified the 
metric around the black hole, the black hole mass and radius that we have obtained coincide with those used by Saida 
and Soda ||^, who assumed a Schwarzschild-like metric. This shows that their ansatz of the Schwarzschild-like metric 
does not introduce a specialization of the problem. 

Let us make a discussion about the definition of mass. In general it is not obvious which definition is the most 
appropriate among many definitions of quasi-local mass. However, in a spherically symmetric and asymptotically flat 
spacetime, the Misner-Sharp mass coincides with the ADM mass at spatial infinity, and with the Bondi-Sachs mass 
at null infinity, as Hayward proved fl^ . Although the asymptotic region in our model is not flat but expanding, 
if the vacuum region inside S is much larger than the black hole, we can regard the regions just inside S as quasi- 
asymtotical-flat and hence the above coincidence makes sense. Among the three the Misner-Sharp mass is the simplest 
to calculate. 

We have not investigated the behavior of the interior vacuum region, which was partially investigated by Scheel et 
al. 1^ and Jacobson Their analyses suggest that the scalar field increases at the horizon and hence the horizon 
area decreases with time as the background universe expands. Although Jacobson's particular solution of the scalar 



5 



wave equation, with very fast variation of (j>{t, oo) oc t in the background universe, look unreahstic, we speculate 
that his conclusion may be unchanged for a more realistic background model because homogenisation of BD field 
occurs more effectively for slower time- variation of (/<(i, oo) [Q. These studies, however, indicate that there are still 
uncertainties to be resolved as to the behaviour of black holes in background universes with arbitrary time variations 
in G. A detailed numerical study of a spherical collapse in a background Friedmann universe will answer some of 
these outstanding questions. 
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